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ABSTRACT 

We present X-ray and optical data on the Be/X-ray binary (BeXRB) pulsar IGR JO 1054- 
7253 = SXP1 1 .5 in the Small Magellanic Cloud (SMC). Rossi X-ray Timing Explorer (RXTE) 
observations of this source in a large X-ray outburst reveal an 11.483 + 0.002s pulse period 
and show both the accretion driven spin-up of the neutron star and the motion of the neutron 
star around the companion through Doppler shifting of the spin period. Model fits to these data 
suggest an orbital period of 36.3 + 0.4d and P of (4.7 + 0.3) x 10" 10 ss _1 . We present an orbital 
solution for this system, making it one of the best described BeXRB systems in the SMC. The 
observed pulse period, spin-up and X-ray luminosity of SXP11.5 in this outburst are found 
to agree with the predictions of neutron star accretion theory. Timing analysis of the long- 
term optical light curve reveals a periodicity of 36.70 + 0.03d, in agreement with the orbital 
period found from the model fit to the X-ray data. Using blue-end spectroscopic observations 
we determine the spectral type of the counterpart to be O9.5-B0 IV- V. This luminosity class 
is supported by the observed V-band magnitude. Using optical and near-infrared photometry 
and spectroscopy, we study the circumstellar environment of the counterpart in the months 
after the X-ray outburst. 

Key words: X-rays: binaries - stars: spectral classification, Be - ephemerides - Magellanic 
Clouds 



1 INTRODUCTION 

Be/X-ray binaries (BeXRBs) are systems in which a compact ob- 
ject orbits a massive, early type star that at some stage has shown 
evidence of line emission in the Balmer series. These stars are main 
sequence stars with a luminosity class of III to V and as such are 
a different sub-class of high-mass X-ray binary (HMXB) to the su- 
pergiant systems. T o date there are no BeXRBs that contain a con- 
firmed black hole ( iBelczvnski & Ziolkowski] l2009h meaning the 
detection of pulsed X-ray emission (caused by the strongly magne- 
tized neutron star; NS) is an extremely robust tool for identifying 
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such systems. The NS accretes via interactions with an extended 
envelope of material in the equatorial plane of the Be star. These 
systems typically have wide, eccentric orbits meaning most X-ray 
outbursts are due to the NS passing through periastron where the 
density of accretable material is greatest. These outbursts are de- 
noted Type I outbursts and are generally in the luminosity range 
10 36 - 10 37 erg s _I and last from a few days to a few weeks, depend- 
ing on the binary period of the system. The less common Type II 
outbursts are brighter, > 10 37 erg s~', and have no correlation with 
orbital phase. The cause of these outbursts is most likely an en- 
larged circumstellar disk that has grown to encompass part or all of 
the NS orbit, meaning accretion o ccurs for prolonged periods; fr om 
several days to a few months (see Stellar, W hite & Rosnerl 19861 for 
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further details of these outburst types). During phases of accretion, 
the pulse period of the pulsar is often seen to decrease, suggesting 
that large transfer of angular momentum from the companion to 
the NS is taking place. This is known as spin-up. This phenomenon 
can be seen on both long-term (several years) and short-term (sev - 
eral days) time scales (e.g. SXP46.6, SXP144; lGaiache et al.l2008T) . 
Long-term changes in spin period are a complex mixture of spin-up 
from epochs of X-ray activity and spin-down from epochs of X-ray 
quiescence. Short term changes are most apparant under larger ac- 
cretion rates and hence, high rates of angular momentum transfer. 
Therefore, observing these systems during a Type II outburst pro- 
vides the best data to study accretion driven spin-up. 

Interactions between the NS and the companion star are 
clearly important in determining the X-ray activity of any particular 
system, and so studying the circumstellar disk around the main se- 
quence star is becoming an ever more prominent field of research. 
Optical monitoring by the Optical Gravitational Lensing Experi- 
ment (OGLE) on the 1.3m Warsaw telescope at Las Campanas ob- 
servatory, Chile, has provided several years of I-band photometry 
of Magellanic Cloud stars. The relatively cool temperatures of the 
circumstellar disk in Be stars means that this material contributes 
strongly to the continuum emission in the near infrared (NIR). 
Thus, variability in the OGLE light curves of BeXRBs is thought 
to indi cate changes in t he stru ctur e or size of the circumste llar ma- 
terial. I^wnsendetd] {2Ql3) and iMcGowan etalj d2008l) present 
examples of how flares in the optical light curves of some BeXRB 
systems can be associated with enhanced X-ray activity. In many 
cases, the variability is shown to be periodic or quasi-periodic; the 
latter thought to be associated with the growth and decay of the 
Be star's circumstellar disk. Temporal analysis often shows the pe- 
riodic variability is caused by the presence of a compact object. 
The orbital periods of several BeXRBs have been discovered, or 
found to agree with a previously derived X-ray period , in this way 
(e.g. ISchurch et ai]|20ld : ISchmidtke & Cowlevl [2006). Variations 
in the size of the disk can a lso be studied through Ha spectroscopy. 
iGrundstrom & Giesl {2006) present numerical models of circum- 
stellar disks and show that the strength of the Ho- equivalent width 
is intrinsically related to the radius of the circumstellar material. 
These models rely on several observable parameters including the 
temperature of the star and the inclination of the disk to the ob- 
server; parameters that are obtainable by knowing the spectral type 
of the star and from the careful analysis of X-ray data like that out- 
lined above. 

The Small Magellanic Cloud (SMC) is home t o nearly 60 
BeXRB systems dCoe et alj|2008l; ICorbet et al.ll2008T) as well a s 
several candidate HMXB systems dShtvkovskiv & Gilfanov 20051) . 
The system that is the subject of this paper, IGR J01054-7253, was 
discovered as a new X-ray source in th e direction of the wing of the 
SMC in June 2009 dBozzo et al]|2 009). It was later confirmed as a 
BeXRB system (see section 2 for details). In this paper we present 
extensive X-ray and optical data of this XRB. Section 2 gives de- 
tails of the X-ray observations made and presents X-ray data of 
the outburst. Section 3 describes the orbital model fitting to sev- 
eral weeks of RXTE data and presents the orbital solution for the 
system. Section 4 presents optical and IR data of the companion 
star. In section 5 we use Science Verification (SV) data from the 
new broadband spectrograph X-shooter on the Very Large Tele- 
scope (VLT) to spectrally classify the optical counterpart. Section 
6 is a discussion of our findings. We review the implications of the 
accretion properties and orbital solution in a broader context and 
discuss what might be happening to the circumstellar disk in this 
system. We end with our conclusions in section 7. 



Table 1. INTEGRAL observations of the SMC and 47 Tuc in which IGR 
J01054-7253 was detected 



Revolution 


MJD IW „ 


Exposure Time (ks) 


Count Rate (ctss 1 ) 


812 


54989.2 


63.3 


0.4927 


813 


54992.2 


53.1 


1.1790 


814 


54995.2 


60.8 


1.1149 


815 


54998.2 


58.8 


0.9149 


816 


55001.2 


62.1 


0.8718 


817 


55004.2 


63.1 


1.31728 


818 


55007.2 


40.2 


1.35568 



2 X-RAY DATA 

2.1 INTEGRAL 

As part of a key programme monitoring campaign of the SMC and 
47 Tuc, INTEGRAL observed the SMC for approximately 90 ks per 
satellite revol ution (~3 days) between 2008 November 1 1 and 2009 
June 25 (see lCoe et al]|201Cl for more details). IGR J01054-7253 
was detected by an I NTEGRAL/IBIS obs ervation on MJD 54989 as 
a new X-ray source dBozzo et al]|2009l) . Table [T] gives the obser- 
vation dates and exposures relating to the detection of this object. 
Unfortunately, these were the last seven observations of the key 
programme and Rev. 812 was the first observation of the SMC for 
nearly 45 days. As such we cannot be certain how long the source 
had been in outburst for before INTEGRAL discovered it; it was in 
the RXTE field of view, but only at a collimator response of approx- 
imately 0.15, making sources harder to detect. Likewise, we only 
detected the source with INTEGRAL up to Rev. 818 which, as dis- 
cussed below, was nearly two months before the source switched 
off. 

Individual pointings (science windows) were processed 
using the INTEGRAL Offline Science Analysis v.7.0 (OS A) 
( Gold wurm etail l2003h and were mosaicked into revolution sky 
maps using the weighted mean of the flux in the 3-10 keV (JEM-X) 
and 15-35 keV (IBIS) energy ranges. Light curves in these energy 
bands were generated on science window (~2000 s) and revolution 
time-scales. Fig. [T] shows the JEM-X (top panel) and IBIS (mid- 
dle panel) revolution light curves of IGR JO 1054-725 3. The JEM- 
X light curve appears to be mostly constant within errors, whereas 
some variation in the hard X-ray flux is apparant from the IBIS light 
curve. By eye it also appears that the hard and soft X-ray flux are 
anticorrelated with one another. However, plotting the hardness of 
the outburst in the ratio ( 15-35)/(3-10) keV against time showed no 
real evidence for a change in the hardness of the X-ray emission. 
Timing analysis on the INTEGRAL light curves was not possible 
due to low count rates. 

2.2 RXTE &■ Swift 

A Swift/XRT follow up observation of the source took place on MJD 
54998.6 aimed at refining the INTEGRAL position and searching 
for pulsed X-ray emission. The observation revealed a precise po- 
sition for this source of RA 1:04:41.41, Dec -72:54:04.6 (2000) 
with a 3.6 arcsec error circle dCoe et alj f2009). This position lies 
within 3.5 arcsec of the V=14.8 star [M2002] SMC 59977 which is 
now thought to be the optical counterpart. Further data on the coun- 
terpart are presented in sectio ns 4 & 5. Timi ng analysis performed 
on the XRT data is reported in lCoe et"aL (2009). Those authors sug- 
gest that the most significant periodicity was 17.49s, although this 
is now known to be incorrect. The correct 1 1 .48s period was not 
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Figure 1. Combined INTEGRAL and RXTE light curve of IGR J01054-7253. The JEM-X flux (top panel) and RXTE pulsed flux (bottom panel) are both 
in the 3-10 keV band, whilst the IBIS flux (middle panel) is in the 15-35 keV band. The two RXTE points marked with a cross rather than a diamond are 
the observations used to extract spectra - see section 2.3 and Fig. [6]for details. The vertical dashed and dot-dashed lines represent the time of Swift and 
XMM-Newton observations respectively. 
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Figure 2. Power spectrum of IGR J01054-7253 showing the 1 1.5s spin pe- 
riod and the third, fourth and fifth harmonics. The second harmonic is not 
seen. SMC X-1 is also detected at higher frequencies. 



foun d in the XRT data until after RXTE had found the correct pe- 
riod dCorbet et alj|2009t IGR J01054-7253 will here-on be refered 



to as SXP11.5). A Swift spectrum was also extracted and is dis- 
cussed in §2.3. RXTE observations began three days after the Swift 
observation, once the source location had been refined, and lasted 
for just over two months. The PC A made a total of 18 exposures in 
Good Xenon mode during this time, meaning a time resolution of 
just under lyus. The average exposure time was ~6ks per observa- 
tion. These data were then binned at 0.01s before being background 
subtracted to produce cleaned light curves in the 3-10 keV band. 
The cleaned light curves were then barycentre corrected and nor- 
malised to the number of PCU's on at each time interval during 
the observation. The RXTE pulsed flux detections are presented in 
the lower panel of Fig. [T] As can be seen, the pulsed flux peaked 
at around MJD55020 and steadily declined until dropping below 
the RXTE sensitivity level (about 1 x 10 35 ergs~' at the SMC). The 
peak detected luminosity of the source was 1 x 10 37 ergs~' in the 
3-10 keV band, using the average of several pulsed fraction mea- 
surements of 0.16. However, we stress that the pulsed fraction was 



only measured using a simple 



max— mm 



method and does not take into 



account the complex shape of the pulse profile; although the aver- 
age value of 0. 16 from the RXTE observations is in agreement with 
the XMM-Newton measurement (see §2.3) and the pulse shape does 
not change by much during the outburst, suggesting that these mea- 
surements are quite robust. Timing analysis was performed on these 
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Figure 3. RXTE 3-10 keV light curves folded at the observed spin period 
during the peak of the outburst (top panel: MJD 55019) and just before 
the end of the outburst (bottom panel: MJD55056). The pulse profiles have 
been smoothed and arbitrarily shifted in phase for viewing purposes. 
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Figure 4. RXTE pulse profiles in the soft (3-10 keV) and hard (10-30 keV) 
energy bands from the observation on MJD 55023. 



data using a Lomb-Scargle periodogram to search for the spin pe- 
riod of the NS. An example power spectrum that is representative of 
most of the observations is shown in Fig. [2] The dominent period- 
icity detected was the third harmonic, with the fundamental period 
and the fourth and fifth harmonics detected at lesser powers. This 
behaviour is also seen in the phase-folded light curves. Fig.[3]shows 
two examples of the 3-10 keV light curves folded at the spin period 
observed by RXTE. The profile in the upper panel shows how the 
third harmonic dominates the emission on MJD 55019 with an un- 
usual triple structure; this shape is representative of all of the other 
profiles until the last few observations of the outburst. The profile 
in the lower panel shows how the emission geometry had changed 
into a profile dominated by the fundamental frequency by the end 
of the outburst, as shown by a much more single peak dominated 
profile. We also note here that a search over the entire 12 years of 
monitoring data from RXTE reveals no previous detections of this 
source, although this may be due to a consistently low collimator 
response of less than 0.2. Another way of exploring the emission 
on the NS surface is to divide an observation up into several en- 
ergy bands and fold each light curve on the detected period, giving 
energy dependent phase-folded light curves. Fig.|4]shows the 3-10 
keV (soft) and 10-30 keV (hard) RXTE folded light curves from 
the observation made on MJD 55023. The soft profile is very simi- 
lar to that shown in the top panel of Fig. [3] demonstrating that there 
was little variation in the emission in this energy range. Comparing 
the soft and hard profiles however, shows evidence for a change in 
structure; the 3rd peak in the profile changes from being smaller 
than the other two to a comparible height. This is possible evidence 
that the emission mechanism is different at different energies. This 
will be discussed further in section 6. 
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2.3 XMM-Newton 

SXP11.5 was sere ndipituously detect ed during one observation of 
the X MM-Newton jjansen et alfeOQlh large programme SMC sur- 
vey dHaberl & Pietsch1l2008r) on MJD 55011.24-55011.59 (ObsID 
060121 2601). The source was visible only in CCD7 of EPIC- 
MOS1 dTurner et alj|200lh at the very rim of the detector, where 
a light path outside the nominal fi eld of view (FoV) e xists due to a 
gap in the instrument structure fcf. lTurner et alj200l[ Fig. 11). We 



Figure 5. XMM-Newton pulse profile of SXP1 1.5 derived from MOS 1 data 
in different energy bands. The profiles are background subtracted and then 
plotted as a ratio to the mean count rate. The mean count rates for the energy 
bands used are (from top to bottom): (1.24, 3.44, 2.88, 1.18 and 8.75) x 10~ 2 
cts/s. 
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note, that the calibration for this region might not be as advanced 
as for the rest of the detector. 

We used XMM-Newton SAS10.0.(Q for data processing. We 
identified sources in the FoV for astrometric boresight correction 
(ARA=-3.49", ADec=0.63") and obtained for SXP11.5 a revised 
XMM-Newton position of RA 01:04:41.60, Dec -72:54:04.5, which 
agrees very well with the position determined by Swift. For the sys- 
tematic error we estimate ~4 arcsec, since the source lies at a large 
off-axis angle of ~17 arcmin. 

For data reduction we selected events with PATTERNS 12 in an 
elliptical extraction region, placed on the source, and a box, lying 
on the point source free part of the light-leak region, respectively. 
Filtering of periods with high background was not necessary, since 
soft proton flares were at a quiescent level, yielding a net deadtime 
corrected exposure of 24.3 ks. For spectral analysis we used only 
events with flag = Sxfffeffff (similar to flag=8, but includ- 
ing events outside the nominal FoV). The spectrum was binned to 
contain ^20 counts bin -1 and response matrices and ancillary files 
were created using the SAS tasks rmfgen and arfgen. 

The power density spectrum of this observation exhibits a 
peak at the period found by RXT E at 0.087 Hz. Using a B ayesian 
periodic signal detection method ^Gregory & Loredol 19961) yields a 
pulse period of (1 1.483 ±0.003) s (ler error). The XMM-Newton de- 
tection independently confirms both the position of this source and 
its spin period. The folded pulse profile in different energy bands 
is plotted in Fig. [5] The shape of the 0.2-10.0 keV profile is rel- 
atively sinusoidal and a calculated pulsed fraction of 0.16 agrees 
with RXTE measurements. Although the third harmonic is not seen 
in the XMM-Newton data, probably because it is close to the de- 
tector timing resolution of 2.6s, the profiles show that there is very 
little variation in the X-ray emission amongst these energy bands. 

Fig. [6] presents the simultaneous fit to the XMM-Newton 
(black) and RXTE spectra (red and green taken 10 days before and 
12 days after the XMM-Newton detection respectively - see Fig.[T](. 
The two RXTE spectra were chosen as they are the ones taken clos- 
est in time to the XMM-Newton spectrum that contain no other ac- 
tive pulsar in the field of view, based on the RXTE power spectrum 
for that observation. The spectra were modelled with an absorbed 
power-law with a high-energy cutoff allowing only a constant nor- 
malisation factor between the three spectra with the assumption of 
no change in the spectral shape (only intensity changes). The flux 
during the XMM-Newton observation was apparently significantly 
lower (although there is some uncertainty due to the far-offaxis 
source position). We fix ed the Galactic photoelectr ic absorption 
at an N H = 6xl0 20 cm- 2 jDickev & Lockmanj [l990). whereas the 
SMC column density with abundances at 0.2 for metals was a 
free fit parameter. We obtained a best-fit with^ 2 /dof = 184/179 
with the best-fit parameters: N h>S mc = 2.30+°^ x 10 21 cnr 2 , T = 
1.06!°;^, cutoff energy = 9.24^ keV, folding energy = 22.0^ 
keV and a detected flux of (7.81 ± 0.57) x 10~ 12 ergs cirT V in the 
3-10.0 keV band. Assuming a distance of 60 kpc, this corresponds 
to an unabsorbed XMM-Newton luminosity of L3_io.o = (3.35 ± 
0.25) x 10 36 ergs s . The luminosity of the two RXTE observations 
found using the normalisation fit are L 3 _ 10 .o = (1.79 ± 0.07) x 10 37 
ergs s and L3_io.o = (2.54 + 0.11) x 10 37 ergs s" 1 for the red 
and green spectra respectively. These agree nicely with the value 
obtained using the pulsed flux light curve and the pulsed fraction. 
The XMM-Newton value is much lower than the values found in the 
RXTE observations. It is unclear how much of this difference is due 
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Figure 6. RXTE & XMM-Newton spectra recorded during the outburst. The 
XMM-Newton spectrum was taken on MJD 55011 and is shown in black. 
The RXTE spectra were taken on MJD 55001 & 55023 (red and green 
respectively). All three spectra were fit simultaneously with an absorbed 
power-law with high-energy cutoff allowing only a constant normalisation 
factor between the three spectra under the assumption that the spectral shape 
does not change. The model fit is presented in the text. 

to source variability and how much is due to the unusually far off- 
axis position of the source on the MOS detector. For comparison, 
this model was also fit to the spectrum from the Swift observation. 
The luminosity was found to be L 3 _ 10 .o = (9.51 + 0.48) x 10 36 ergs 
s" 1 . This value is a factor of 2 smaller than the RXTE value mea- 
sured just 3 days later, but could be explained by variability within 
the source at the start of the outburst. 

SXP11.5 was also in the FoV of a previous XMM-Newton 
observation on MJD 54010.99-54011.26 (ObsID 0402000101), 
when no significant X-ray emission at this position was detected. 
Analysing the EPIC-pn data, we found a 3cr upper limit of 0.002 
cts s -1 for the 0.2-12.0 keV energy band, which, assuming the same 
spectrum as above, corresponds to an unabsorbed luminosity limit 
of Lo.2-10.0 < 6.4 x 10 33 ergs s _I . 



3 ORBITAL SOLUTION TO IGR J01054-7253 

It became apparent after the first few RXTE observations that there 
were both spin-up and spin-down trends in the pulse period of 
SXP11.5, indicating that we were seeing the motion of the NS 
around the companion star through Doppler shifting of th e spin pe- 
riod, simil ar to that seen in SX P18.3 dSchurch et alj2009h and GRO 
J1750-27 dShawetal .1120091) . Fig. [7] shows the detections of the 
third harmonic of the pulse period during the outburst and the as- 
sociated errors. Overplotted is the model fit to the data as discussed 
later. The third harmonic was used to fit the model because it was 
several times more powerful in the power spectrum than the funda- 
mental frequency (Fig. [2)1 and so the associated errors are smaller. 
The error bars increase in size towards the end of the outburst as 
the pulsed emission became much weaker as the pulsar switched 
off. The final two data points are one third of the detected funda- 
mental period as the third harmonic had disappeared by this point 
(see discussion) and the first data point is the Swift detection, which 
is also one third of the detected period as, like the MOS-1 detector 
on XMM-Newton, we were unable to detect the third harmonic due 
to longer binning needed with such a low count rate source. Due to 
the length of the gaps between observations, it was not possible to 
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Table 2. The orbital parameters for IGR J01054-7253 from the analysis of 
3-10 keV RXTE PCA data. 
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Figure 7. Period of the third harmonic of SXP11.5 as detected by RXTE 
during the two months of pointed observations. The final two data points 
are one third of the detected fundamental period as the third harmonic had 
disappeared by this point and the first data point is one third of the period 
detected by Swift (see text). Overplotted is the model fit to the data as de- 
scribed in the text. 



match phase between them and do more accurate pulse arrival time 
analysis. 

As mentioned, the oscillatory nature of the spin period is most 
easily explained by the Doppler shifting of the X-ray light from the 
NS by the binary orbital motion in the system. However, a general 
spin-up trend, most likely caused by angular momentum transfer 
during the accretion of matter onto the NS, is also apparent in the 
data. Consequently the observed data were fit with a simple spin-up 
model given by 



P(t) = P(t ) + P(t - t ) - 



Ht - t ) 2 



(1) 



where P is the spin period of the neutron star, P and P are the 
spin-up and change in spin-up, and t and t are the start time and 
time since the start time over which the model is fit. This model 
was convolved with a standard orbital model which calculates the 
line of sight velocity of the NS. This combined model iterates over 
some input parameters until a best fit for a Levenberg-Marquardt 
least-squares fit is achieved. Standard values for other HMXB sys- 
tems were used for initial parameters that were unknown. After first 
running a simplified version of this model that assumed a circular 
orbit, it was clear that an acceptable fit could only be achieved by 
including an eccentricity component and excluding P as the data 
were not sufficient to fit the second derivative. The final fit gave 
an orbital period of 36 .3 ± 0.4d and a P of (4.7 ± 0.3) xlO^^ss" 1 
dTownsend e7ai1l2009l) . The model fit to the data is shown in Fig. 
[7] and the various parameter values are displayed in Table [2] along 
with the reduced chi-squared value. 



4 OPTICAL AND IR DATA 

In this section, we present spectroscopic and photometric data in 
the optical and NIR wavebands from a variety of ground-based 
telescopes. The aims of this were to spectrally classify the opti- 
cal counterpart and to explore the changes in the circumstellar disk 
leading up to, during and after the X-ray outburst. 



Parameter 




Orbital Solution 


Orbital period 


P orbital (d) 


36.3 ± 0.4 


Projected semimajor axis 


a x sini (light-s) 


167 ±7 


Longitude of periastron 




224 ± 10 


Eccentricity 


e 


0.28 ± 0.03 


Orbital epoch 


Tperiastron (MJD) 


55034.3 ± 1.0 


Spin period 


P(s) 


11.48143 ±0.00001 


First derivative of P 


P(ss-') 


(~4.67±0.3ir 10 


Goodness of fit 




0.82 



4.1 Observatories and Instrumentation 

• OGLE III - The OGLE project has been steadily monitoring 
millions of stars in the Magellanic Clouds f or the past 1 2 years (see 
lUdalski. Kubiak & Szvmafiskil d 19971) and ISzvmanskil d2005l) for 
more details on the OGLE instrumentation and catalogue). Photo- 
metric images in the I-band are taken almost every night whilst the 
SMC is visible; these images now make up a decade long database 
that includes most of the HMXB systems in the SMC. 

• Faulkes Telescope - FT South, located at Siding Spring, Aus- 
tralia is a 2m, fully autonomous, robotic Cassegrai n-type reflec- 
tor w hich employs a Robotic Control System (RCS) lTsapr as"et al.l 
2009). The telescope was used both in Real Time Interface (RTI) 
mode and Offline mode for the observations of SXP11.5. All the 
observations were pipeline-processed which does the flat-fielding 
and de-biasing of the images. 

• 1.9m Radcliffe telescope at the South African Astronomical 
Observatory (S AAO) - Photometry was done using the SAAO CCD 
at the Cassegrain focus. The spectra were obtained using the unit 
spectrograph combined with a 1200 1/mm grating and the SITe de- 
tector at the Cassegrain focus. Data reduction was performed using 
standard IRAF packages. 

• 1.4m Infrared Survey Facility (IRSF) telescope at SAAO - The 
IRSF is a Japanese built telescope designed specifically to take 
simultaneous photometric data in the J, H & K s bands with the 
SIRIUS (Simult aneous three-colour Inf raRed Imager for Unbiased 
Survey) camera (Na gashima et alj fl999). Data reduction was per- 
formed using the dedicated SIRIUS pipeline based on the National 
Optical Astronomy Observatory's (NOAO) IRAF software pack- 
age. The pipeline was provided by Yasushi Nakajima at Nagoya 
University, Japan. This performs the necessary dark subtraction, 
flat fielding, sky subtraction and recombines the dithered images. 
Photometry on the reduced images was performed using standard 
IRAF routines. 

• X-shooter - The first of the secon d generation VLT instru- 
ments, X-shooter (D'Odorico et al . 2006), is a three arm, single ob- 
ject echelle spectrograph for the Cassegrain focus of one of the VLT 
UT's. The instrument simultaneously covers the wavelength range 
300-2400 nm at resolving powers R = ^ = 5100, 8800 and 5600 in 
the UVB (AA = 300-550 nm), VIS (AA = 550-1015 nm), and NIR 
arms (AA = 1025-2400 nm) respectively. However, the resolution 
obtained is dependent on the slit width and seeing conditions. 



4.2 Photometry 

Fig. [8] shows the full 12yr optical light curve of SXP11.5 from 
OGLE monitoring (open black triangles) and coverage from 
Faulkes Telescope (FT; closed red squares) south and the 1.9m, 
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5.3x10* 5.4x10* 5.5x10* 

Date (MJD) 

Figure 8. Combined OGLE III, Faulkes Telescope (FT) south and SAAO 
1.9m I-band light curve of the optical counterpart in SXP1 1.5 (open black 
triangles, closed red squares and open blue square respectively). The FT I- 
band data have been transformed into the OGLE I-band for direct compari- 
son of the magnitude of the star. The vertical dashed lines indicate the start 
and end of the Type II X-ray outburst as seen by INTEGRAL and RXTE; this 
time period spans approximately two orbital cycles, illustrating the scale of 
the optical variability. 




SAAO telescope (open blue square). The FT I-band data were cal- 
ibrated onto the OGLE I-band using 18 field stars from the OGLE 
database that show no variability over the set of measurements. The 
OGLE data are sparse compared to the normal coverage of the SMC 
because this system occasionally falls onto a gap in the detector 
chip which is dependent on its orientation angle. Despite this, the 
data presented were detrended by subtraction of a 50 day moving 
mean model to remove large amplitude variability, allowing tim- 
ing analysis to be performed. A Lomb-Scargle periodogram of the 
detrended light curve shows the presence of a 36.70 ± 0.03d peri- 
odicity above the 99.9% significance level, as shown in Fig. [9] This 
periodicity was only detected in certain sections of the light curve, 
but confirms the orbital period predicted from the X-ray data. Fig. 
[9] also shows the phase-folded light curve folded at the detected 
period and phase alligned to the ephemeris given in Table [2] The 
profile seems to be quite sinusoidal, being brightest just after pe- 
riastron and faintest just after apastron, suggestive of a slight lag 
between the periodic brightening of the disk with the phase of the 
NS. The narrow peak and trough at phases of approximately 0.4 
and 0.9 respectively are somewhat more peculiar; they may suggest 
something is happening half way between the optical minimum and 
maximum. However, on folding the light curve at various phases or 
with different bin sizes these features come and go, suggesting that 
they may not be a real physical feature of the system. There are 
no MACHO or OGLE II data of this source. The OGLE, FT and 
SAAO data were taken before, during and after the X-ray outburst 
respectively (see discussion section for interpretation). The orbital 
period found here was also used as a frozen parameter in a second 
orbital model fit to the X-ray data. This did improve slightly the 



Figure 9. Lomb-Scargle periodogram of the detrended section of the OGLE 
III light curve between MID 53400 and MJD 54200 (top). A peak at 36.70 
± 0.03d is seen above the 99.9% significance level, verifying the orbital 
period predicted by the X-ray data. This period was only detected in this 
section of the light curve when the disk was seemingly at its lowest flux 
level. Analysis of other individual sections and the light curve as a whole 
did not reveal any significant periods. The phase-folded light curve (bottom) 
is folded such that a phase of 1 .0 in the figure corresponds to periastron of 
the NS based on the result of our X-ray fit. 

errors on the parameters presented in Table [2] but only on the order 
of a few percent. 

Photometric data taken on the IRSF and the SAAO 1.9m tele- 
scope are shown in Table [3] along with measurements previously 
published in catalogues. The IR SF Magellanic Clouds Point Source 
catalogue bv lKatoetal.1 i l2007l) was used to calibrate the IRSF data. 
The Magellan i c Clou ds Photometric Survey (MCPS) catalogue by 
IZaritskv et al] (2002) was used to calibrate the 1.9m optical data. 
The optical photometry is not as good as the IR as can be seen 
from the associated errors. This is due to the poor seeing at optical 
wavelengths during the observations, and variability in the condi- 
tions during the long integration times required because of a full 
moon. Also, the smaller field of view of the SAAO CCD compared 
with the SIRIUS CCD meant there were fewer stars which could be 
used to perform the calibration. A pipeline provided by SAAO was 
used to calculate the PSF magnitude of the optical images after per- 
forming flat fielding. The error is calculated by taking the average 
of the square of the magnitude errors from the MCPS catalogue, for 
the stars used for calibration, and the standard deviation of the stars 
from the fitted line. The final error is , la 2 . , + cr?.. Data from 

\ avg.catalogue fit 
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Table 3. IR and optical data of the counterpart to SXP1 1.5. 



Catalogue 


Date (MJD) 


J 


H 




2MASS 1 
SIRIUS 2 


51035 
52894 


14.16 ±0.03 
14.21 ± 0.08 


13.70 ± 0.04 
13.75 ±0.06 


13.52 ±0.04 
13.74 ±0.04 


Telescope 


Date (MJD) 


J 


H 


K s 


IRSF 


55196 


14.23 ± 0.04 


13.75 ± 0.03 


13.61 ±0.04 


Telescope 


Date (MJD) 


B 


V 


I 



1.9m 55200 14.91 ±0.31 14.91 ±0.11 14.59 ±0.10 



ISkrutskie etld] 120061) . fcto et alj fe007l) . 



the IRSF Magellanic Clouds Po int Source catalogue and 2MASS 
catalogue dSkrutskie et alj|2006h are presented as reference to past 
IR observations of this source. Unfortunately, the IR data are not 
directly comparable to the I-band data in the OGLE light curve as 
both measurements fall into gaps in the data and as such we can- 
not say if the NIR bands follow what is happening in the I-band. 
However, we can make two cautious observations; firstly, that there 
is little variation in the J, H & K band values presented in Table [3] 
and, secondly, that at the time of the SIRIUS catalogue measure- 
ment the I-band seems to be much fainter than at the time of the 
IRSF measurement (assuming the I-band at MJD 53000 is around 
15.2 magnitudes in comparison to the present value of 14.6). This 
suggests that there may be a change in the (I-K) colour in this sys- 
tem. If true then this is indicative of a changing temperature within 
the optically thick circumstellar disk. 



4.3 Spectroscopy 

SXP11.5 was observed during the X-shooter Science Verification 
(SV1) phase in August 2009 (MJD 55055; ~10d before the X-ray 
outburst ended), with a total integration time of 1200^. The spec- 
tra were reduced wit h a beta version of the ESO X-shooter pipeline 
( Gold onT et al J20061) . which uses subtract i on of the sky lines based 
on the procedures developed by Kelson (2003). The pipeline re- 
duction used calibration spectra taken during the commissioning 
run for order location and tracing, flat fielding and wavelength cal- 
ibration. The final product from the pipeline is an extracted 2D, 
wavelength calibrated, rectified spectrum with orders combined us- 
ing a weighting scheme. For further data processing and analysis 
we used a combination of IDL programs and standard IRAF tools. 
The spectra were then smoothed with a boxcar average of 7, nor- 
malised to re move the continuum and then shifted by 150 km s 
dAllenll 19731) to account roughly for the recession velocity of the 
SMC and hence to place spectral features at approximately the cor- 
rect wavelengths. These spectra are presented in Figs. 1 101 and 1 121 
The UVB arm spectrum is used in the next section to spectrally 
classify SXP1 1.5, whilst the VIS arm spectrum is used to measure 
the Ha equivalent width just before the outburst came to an end. 
This is discussed and compared to a SAAO Ha measurement taken 
months later in section 6. 



2 Program ID 60.A-9439(A) 



5 SPECTRAL CLASSIFICATION 

Classification of Galacti c Be stars relies on using t he ratio of 
many metal-helium lines ( Walborn & Fitzpatrick 1990). However, 
this type of classification, based on the Morgan-Keenan (MK; 
Morg an. Keenan & Kellmanll 19431) system, is particularly difficult 
in the low metallicity environment of the SMC because these metal 
lines are either very weak or not present a t all. Usi ng high signal- 
to-noise ratio spectra of SMC supergiants. lLennor] J 19971) devised 
a system for the classification of stars in the SMC that overcomes 
the problems with low metalicity environments. This system is 
normalized to the MK system such that stars in both systems ex- 
hibit the same trends in their line strengths. Another difficulty is 
that the Balmer lines in particular will be rotationally broadened 
due to the high rotational velocities of Be stars and hence may 
obscure any comparisons to closely neighbouring lines. We have 
thus used the classification method as laid out in iLennonl d 19971) 
and utilised further in Eva ns et all d2004l) . For the luminosity clas- 
sification we have adopted the classification method set out in 
IWalborn & Fitzpatrick! ( 1199(1 . 

The classification of SXP11.5 has been made using the X- 
shooter UVB arm spectrum presented in Fig. [TO] Immediately ob- 
vious is the presence of ionised Helium in the spectrum: He II 
/M4686, 4 541 are present, meaning the sp ectral type must be B0 
or earlier dLennonlfl997l; lEvans et alj|2004l) . He II i4200 is weak, 
and much weaker than He I /14143, meaning it is later in type than 
09. There may be evidence of some Si lines in the spectrum which 
could stretch the classification to Bl, although these are very diffi- 
cult to confirm above the general SNR of the spectrum. Given the 
presence of He II /14686 it would seem that, even if the metal lines 
are supressed by the low metallicity of the SMC or masked by ro- 
tational broadening of t he Balmer and helium lines, a spectral type 
as late as B 1 is unlikely. IWalborn & Fitzpatrick! ( fl990h also present 
spectra with apparent Si lines and classify them as B0, indicating 
that a classification of B0 would be correct in this case if the Si lines 
are present. Given that the metal lines in the spectrum are not so ob- 
vious, we were limited to using the He I /14121/He I /I4143 r atio to 
determine the luminosity class dWalborn & Fitzpatrick!! 199d) . This 
line ratio strengthens towards more luminous stars, suggesting that 
this star has a luminosity class of IV- V. We have also performed 
a check on this luminosity classification by comparing the abso- 
lute magnitude of the source in the V-band with a distanc e modulus 
for the SMC of 18.9 dHarries, Hilditch & Howartr]|2003l) , to deter- 
mine whether the absolute magnitude of the source is consistent 
with the estimated luminosity class for t his spectral typ e (using ab- 
solute magnitudes for OeBe stars from Wegner 2006). A V-band 
magnitude of 14.9 (Table [3} confirms a luminosity class of IV-V 
for a spectral type of O9.5-B0, as a star of higher luminosity class 
would need to be of a later spectr al type. Our spectra l classification 
agrees with recent work done by Mas etti et alJd20l5T) on the SAAO 
1.9m telescope, although those authors give a luminosity class of 
III which may be too early based on the available line ratios in the 
spectrum. 



6 DISCUSSION 
6.1 X-ray behaviour 

Strong triple peaked structure is visible in the X-ray pulse profile 
when the 3rd harmonic was strongest (top panel Fig.[3}- This triple 
structure remained constant for several observations whilst the out- 
burst was at its most intense. Only once the flux had significantly 
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Figure 10. Spectrum of SXP11.5 taken with the UVB arm of the X-shooter spectrograph in the wavelength range 300-550nm. The spectrum has been 
normalised to remove the continuum and the data have been redshift corrected by -150 km s to account for the recession of the SMC. Overplotted are various 
atomic transitions that are significant in the spectral classification of an early type star at their rest wavelengths; He II, He I and metal transitions are in green, 
red and blue respectively. 



dropped did the pulse shape begin to evolve into a more familiar 
single peaked structure (bottom panel Fig. [3}. At this point the 3rd 
harmonic fades and the fundamental begins to dominate the pulsed 
emission. The energy dependent pulse profiles plotted in Figs.|4]and 
[5] demonstrate a possible change in the X-ray emission at higher 
energies. The XMM-Newton profile seems to show the shape of the 
folded light curves are very similar within the 0.2-10 keV energy 
range (although we are aware that some information is being lost 
here due to the none-detection of the 3rd harmonic). This was also 
the case for the RXTE observations in which the 3-10 keV light 
curve was split into smaller energy bands (these have not been pre- 
sented here). However, moving to higher energies in the RXTE pro- 
files we see a change occuring around 10 keV. Below this energy, 
the pulse profile suggests a combination of a pencil beam from one 
pole and a fan beam from the other pole, with an orientation that 
places two sides of the fan beam and the pencil beam approximately 
one third of a phase apart. The smaller of the three peaks is likely 
the pencil beam in this picture, being a phase of 0.5 apart from 
the centre of the fan beam 'double' peak. The hard energy profile 
shows evidence for a number of possibilities; the pencil beam has 
increased in strength, becoming comparible to the fan beam. The 
pencil beam may be narrower, although it remains 0.5 phase away 



from the centre of the other 2 peaks. The fan beam may have a 
narrower opening angle at higher energies given the apparant in- 
filling of the gap between the two peaks. At the moment, these are 
all speculations based on observation and thus more detailed mod- 
elling is required to take this discussion further. 

The orbital para meters in Table [2] are typical o f those seen in 
other BeXRB systems ( Ok azaki & Negueruelal200lh . Most Be sys- 
tems haveP orf , ~ 10-100 d and orbital eccentricities of between 0.3- 
0.5 jBildstenetal]|l997l) . These parameters also pl ace SXP11.5 
firmly in the BeXRB regime of the Corbet diagram dCorbet et al.l 
1986), lending support to the optical determination that this source 
is a BeXRB and not a supergiant system. 

lGhosh&LamblJl979l) show that there is a relationship for bi- 
nary pulsars such that -P a PLi for a given NS mass and magnetic 
moment, where P is the spin-up in sec/yr, P is the spin period in sec- 
onds and L is the X-ray luminosity in units of 10 37 ergs/s. This rela- 
tionship has since been verified bv lCoe. McBride & Corbet! d2010h 
using a much larger database of SMC binary pulsars. The spin pe- 
riod and P from our fit were used along with the highest RXTE 
luminosity t o place this new system onto th e plot of SMC pulsars 
presented in ICoe. McBride & Corbetl (2010). The result is shown 
in Fig. [TT] The position of this source on the plot seems to follow 
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Figure 11. The distribution of SMC BeXRB pulsars on the PvsPLi plane 
adapted from Coe, McBride & Corbel] l2010l) . where P is the spin-up in 
sec/yr, P is the spin period in seconds and L is the X-ray luminosity in units 
of 10 37 ergs/s. SXP1 1.5 is shown to sit in the low period end of the distribu- 
tion, confirming that the luminosity, pulse period and spin- up calculated in 
this an alysis agree with the equations orginally presented in iGhosh & Lamtl 
1 1979) and that this source is, in this sense, normal with respect to the gen- 
eral distribution of SMC pulsars. 



the general distribution of SMC pulsars and thus agrees with the 
theoretical laws of binary accretion. 



6.2 Optical behaviour 

Due to the lower temperatures believed to be present in the circum- 
stellar environment, the photmetric I-band data presented in Fig. [8] 
are a direct measure of the circumstellar emission. These data show 
that this system has a highly variable circumstellar disk and that the 
X-ray outburst was triggered by the disk being in an exceptionally 
large state, relative to its base level. The vertical dashed lines in 
Fig. [D show the start and end points of the X-ray outburst in rela- 
tion to the optical activity. A single observation at SAAO in Jan- 
uary 2010 (MJD 55200) shows that the disk has started to decrease 
in size, although the X-ray outburst had finished before the disk 
had returned to the level at which the X-ray emission began. This 
could be evidence of a lag between the X-ray and optical emission 
in this system. Another puzzle that these data introduce is why we 
only see the orbital modulation in the light curve at the very lowest 
flux level. This would suggest that the variability is coming from 
the stellar a tmosphere as the dis k would seem to be at a minimum. 
In contrast lSchurch et al.l J20101) only detect the orbital modulation 
from SXP18.3 at the very brightest part of the detrended optical 
light curve, suggesting that in this system, the variability is coming 
from the circumstellar disk and not the stellar atmosphere. 

It is currently believed that the amount of red and IR emis- 
sion from the circumstellar disk is directly related to the strength of 
emission lines in this regime. Fig.[T2]shows two spectra of the Ha 
region taken during and many months after the X-ray outburst. The 
top panel is a VLT/X-shooter VIS arm spectrum showing clearly 



Ho- being in emission. This corresponds to the peak of the I-band 
light curve shown in Fig. [8] Equivalent width measurements of 
the line profile give a value of -8.0±0.5A. The bottom panel is a 
spectrum taken on the SAAO, 1.9m telescope 4 months after the 
outburst had ended. This observation was near-simultaneous to the 
SAAO I-band measurement plotted in Fig. [8] The poor SNR in this 
spectrum made it difficult to make a precise equivalent width mea- 
surement, the best estimate being -3 ± 2, so we cannot say by how 
much the disk has reduced in size. However, we can conclude that 
both the Ha and I-band measurements suggest the disk is shrinking. 
Comparision to Ho- measurements of other BeXRBs shows that, 
even during a long Type-II X-ray outburst, the disk in this system 
was comparably quite small (we have equivalen t width measure - 
ments of more than -55 A in some other systems: Icoe et al]|2005h . 
although it s eems not abnormally small based on the relationship 
presented in I Antoniou et alj J2009I) ; those authors present the Ha 
equivalent width of a sample of BeXRBs against their orbital pe- 
riod and show a linear trend. Although SXP 11.5 lies below the best 
fit line of this trend, it is not a complete outlier and may just lie 
within the distribution. Unfortunately a more quantitive handle on 
the IR excess seen in the photometry is needed to help verify the 
true size of the disk relative to other systems in the SMC. So given 
only a small circumstellar disk to accrete from, how does accretion 
continue for over two months in this system? The most plausible 
explanation seems to be that an accretion disk has formed around 
the NS, allowing accretion to continue after the NS has moved away 
from the circumstellar disk which may not, in this case, be filling, 
or be near to filling the orbit of the NS. However, it is not easy 
to detect the presence of an accretion disk in HMXBs and so an 
unusual orbit or other such geometries cannot be ruled out as an 
explanation of what is seen, although this seems unlikely based on 
the fairly standard parameters obtained in the orbital fit of the X- 
ray data. It is also worth noting that the narrow, single peaked Ha 
and Hy6 emission lines seen in the spectra presente d here sugges t 
the disk has a low inclination angle to the observer dStruvdll930h . 
Comparing this to the value of a v sim presented in Table [2] suggests 
that the semimajor axis of the disk could be over 300 light-s (for i 
< 30°), approximately 18 R* assuming a B0 star of radius equal to 
7 R . 

In order to explore the extent of the emission from the circum- 
stellar disk, the contemporaneous optical & IR photometric data 
from MJD 55196-55200 shown in Table [3] were used to compare 
with the predicted stellar emission. The photometric data were first 
dereddened by the establi shed value to the SMC of E(B-V)=0.08 
jSchweriruT & Israel 1 99ll). These were then compared to a Kurucz 
model atmosphere jKurucz|[l979l) for a B0V star (r e// =30,000K 
and log g=4.0). To make the comparison, the model atmosphere 
was normalised to the lowest I-band point in the OGLE III data 
(taken around MJD 54000) on the assumption that this represents 
the occasion on which any contribution from the circumstellar disk 
is at its lowest. The results are shown in Fig.[T3]from which it is im- 
mediately apparent that at the time of the optical & IR photometry 
there was a clear disk contribution across the whole of the optical & 
IR regime. Even the B & V points clearly contain some element of 
such a component. In trying to quantify this excess, we normalised 
the Kuruscz model to the distance estimate of the SMC (60 kpc) and 
found that it lay very close to the I-band normalised model. Work- 
ing from this we converted the I-band flux at this base level into 
a distance estimate, assuming that all the flux is coming from the 
star. We find a distance to SXP 11.5 of 71 ±5 kpc. The error comes 
from calculating the maximum and minimum distance based on a 
range of E(B-V) values. Of course, this value is very dependant on 
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Figure 12. Spectra of SXP11.5 taken in the wavelength range of Ha at 
the VLT (top) and SAAO,1.9m (bottom). The VLT spectrum was taken 
on 2009-08-12 (MJD 55055), just before SXP11.5 fell below the detec- 
tion threshold of RXTE, therefore showing the disk during an X-ray active 
phase. The SAAO spectrum was taken on 2009-12-12 (MJD 55177), four 
months after the X-ray outburst ended. 

obscuration localised to the SMC itself and should be treated with 
some caution. What this calculation shows is that it is very likely 
that the I-band measurement used in the normalisation represents 
very closely the base level of the disk emission in this system. A 
Hff spectrum is needed to confirm this prediction should the I-band 
magnitude ever return to this level. 
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Figure 13. The optical & IR photometry (open squares) from MJD 55196- 
55200 dereddened and compared to a Kuruscz model atmosphere for a B0V 
star (continuous line). The filled triangle indicates the lowest OGLE III I- 
band point which has been used to normalise the Kuruscz model - see text 
for details. This demonstrates the extent of the continuum emission from 
the circumstellar disk in this system. 



7 CONCLUSIONS 

We have presented the complete orbital solution to a newly dis- 
covered high-mass binary system in the SMC, making it one of 
the few well described systems known outside the Galaxy. These 
parameters are similar to other binary systems in the Galaxy and 
Magellanic Clouds and place SXP1 1.5 firmly in the BeXRB region 
of the Corbet diagram. Under huge accretion torques, it is shown 
that the NS is being spun up and that this satisfies current theoreti- 
cal predictions of the relationship between spin-up, spin period and 
luminosity in binary systems. However, it is clear that there is still 
much to learn about the accretion process itself as shown by the 
complex triple peaked structure and evolution of the pulse profiles, 
which are still poorly understood. 

Spectral analysis of the companion has allowed a classifica- 
tion of O9.5-B0 IV-V to be made. Timing analysis of the optical 
light curve confirms the orbital period detected in the X-ray. Con- 
temporaneous Ha spectra show, along with the optical light curve, 
that the circumstellar disk is shrinking in the time since the X-ray 
outburst. Although there is evidence that the disk is unusually small 
compared to other BeXRBs, making the extended period of accre- 
tion seen in this outburst difficult to explain. We show that the low- 
est flux seen in the optical light curve is likely to be from a period 
in which the Be star has entirely lost its disk and make a distance 
estimate to the source based on this assumption. Ha spectroscopy 
will be necessary to confirm this prediction should the I-band flux 
ever return to this level. 

SXP11.5 has proven to be a very interesting addition to the 
SMC binary population; on one hand, it seems to be a very normal 



system showing orbital parameters and a spectral type that match 
very well with previously studied systems. However, it also shows 
peculiarities, such as showing a small Ha equivalent width relative 
to other systems with similar orbital periods. Should this system 
go into another X-ray outburst in the future, further simultaneous 
optical and X-ray measurements are essential to uncovering where 
this system lies in the overall population of SMC BeXRB systems. 
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